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We investigate theoretically acoustic phonon induced decoherence in quantum dots. We calculate the dephas- 
ing of fundamental (interband or intraband) optical transitions due to real and virtual transitions with higher 
energy levels. Up to two acoustic phonon processes (absorption and/or emission) are taken into account simul- 
taneously in a non-perturbative manner. An analytic expression of acoustic phonon induced broadening is given 
as a function of the electron-phonon matrix elements and is physically interpreted. The theory is applied to the 
dephasing of intersublevel transitions in self-assembled quantum dots. 



Understanding decoherence in semiconductor quantum 
dots (QDs) is of prime importance in the context of their po- 
tential application in quantum information. Dephasing of the 
interband transitions in semiconductor QDs have been exten- 
sively studied in the past few years both experimentally and 
theoretically. It has been show that the interaction with lattice 
vibrations (phonons) is a major source of dephasing^l2Jof the 
optical polarization. From the theoretical point of view, inter- 
action between the electronic levels and acoustic phonon can 
be decomposed into diagonal and off-diagonal terms. The di- 
agonal part gives rise to acoustic phonon sidebands^. This 
sidebands are responsible of a rapid and partial decay of 
the polarization which is now well understood-^a. Differ- 
ent mechanisms have been put forward in o rder to explain 
the broadening of the zero phonon line (ZPLj^MUl. I n par- 
ticular, Muljarov and Zimmermann 7 have proposed a mecha- 
nism of virtual transitions triggered by the off-diagonal acous- 
tic phonon interaction. 

More recently, dephasing of the s-p intersublevel transition 
has been studied in QDs^-4 While at low temperature (10K) 
dephasing has been shown to be limited by anharmonic po- 
laron decay mechanism 12 , at higher temperature both real and 
virtual transitions between the two p states were identified as 
responsible for the broadening of the ZPL. 

Virtual transitions refer to two-phonon processes of simul- 
taneous absorption and emission of phonons with same energy 
but different momentum (i.e. elastic diffusion of phonons) 
due to non-resonant phonon couplings (off-diagonal electron- 
phonon interaction) with an higher excited state. This process 
contributes to the decay of the coherence in addition to other 
dephasing processes such as real transitions. 

In self-assembled quantum dot, phonon absorption or emis- 
sion is efficient only if the energy separation between levels is 
of the order of a few meV, corresponding to acoustic phonon 
wavelength of the order of the dot size. However, energy sepa- 
ration between electronic levels is generally larger so that real 
transitions (single phonon absorption or emission) are in gen- 
eral not efficient. This effect is known as the acoustic phonon 
bottleneck effect 13 . Virtual transitions instead (because of 
their non-resonant character) are less sensitive to the energy 
detuning between electronic levels and provide a contribution 
to the dephasing in a larger range of electronic energy separa- 
tion. 

To date, two-phonon processes in QDs have only be treated 
perturbatively with respect to the electron-phonon matrix ele- 



ments. Takagahara has used a self-energy expansion 2 , while 
Muljarov and co-workers have used a cumulant expansiorPEi] 
In the present article, we present a calculation of the dephas- 
ing in quantum dots taking into account up to two-phonon 
processes in a non-pertubative manner. Contrary to previ- 
ous approach, our calculation includes a resummation of all 
the diagrams that involve states which differs from the initial 
reservoir by 1 or 2 p honons. This is achieved using projec- 
tion operator j 15 l 16 l 17 l twice consecutively. An analytic formula 
is obtained as a function of the electron-phonon matrix ele- 
ments, which is physically interpreted and allows to identify 
the contribution of the different phonons to the dephasing pro- 
cesses as a function of their energy. 

I. THEORY OF TWO-PHONON PROCESSES INDUCED 
DEPHASING 

We consider a three level system: the ground state is de- 
noted \g), while |e ) and \ei) are the first and second excited 
states. We study here the coherence of the transition between 
\g) and |eo). We consider the following Hamiltonian: 

H = H e + H ph + Ve-ph (1) 

where H e is the electronic part, H v h is the acoustic phonon 
Hamiltonian, and V e . P h represents the electron-phonon inter- 
action. 

H e = e g \g)(g\ + £ |e )(e | + (s + A)|ei)(ei| (2a) 

Hph = ^ £ i a q a q ( 2b ) 



n.m q 

where the sums over n and m run over g, eo and e\. In quan- 
tum dots, the transition energy e — e g 1S usually large com- 
pared to the energy separation A between the first two excited 
states |erj) and \e\) (in both intraband and interband cases). As 
a consequence, the off-diagonal acoustic phonon couplings af- 
fecting the ground state (Af| e ° and M| ei ) can be neglected. 
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In addition, diagonal interaction in the g state (e. g. for in- 
traband transitions) is removed through a Huangh-Rhys trans- 
formation of the form a q — > a q — M 99 * /e q . The new phonon 
modes are thus the eigenmodes of the crystal vibrations when 
the system is in the ground state \g). In the following, the to- 
tal Hamiltonian will be split as H = Hq + V^_ p h where the 
non-interacting part is Hq = H e + H p ^. 

We are interested in the decay of the optical polarization of 
the g-eo transition due to the interaction of the phonon reser- 
voir assumed to be at thermal equilibrium. We consider an 
initial eigenstate of the form \g,<fro), where \<po) denotes an 
initial phonon state |0o) = \n>i,ri2, ■■■,n>N cr ) where is the 
occupancy of the LA-phonon mode i (rii — ((f>o\ a t &i|<^o)) 
and N cr is the number of atom in the crystal. The limit of an 
infinite crystal is considered {N cr — > +00). In this limit the 
integer values of the n, are distributed according to the Bose- 
Einstein function N(e) = \j(e e l kBT — 1). The dipole oper- 
ator of the g-eo transition is defined as d = \g) (eo| + |eo) (g\. 
The linear polarization (the linear component of the polariza- 
tion in response of a delta pulse excitation) is given bjEJ 



P(t)=i(g,<p \d(t)d\g,<j>o)) 



(3) 



where d(t) = e im / n de- iHt / n . 

As |eo, 4>o) is an eigenstate of the full Hamiltonian H, Eq.[3] 
reduces to: 



P(t)=ie fe «*/»( eOi o | e -iH*/»| eOj o ) 



(4) 



It is known that diagonal acoustic phonon interaction in 
QDs leads to phonon sidebands that are responsible for a 
partial polarisation decay. A Huangh-Rhys transformation 
is made with new phonon modes defined by b q — a q + 
M q ° e °* je q . The new occupancy states are then denoted 



by: 



(&^)"*/vn!|0). As shown in Appendix, P(t) is given 



P(t) = ie-™g(t)P z (t) 



Pz{t) 



(e , 



b \e- im l h \e , 



(5) 



(6) 



where git") is defined in the appendix and corresponds to 
the polarization decay due to the sidebands only, |0o) = 
\n[,n2, nN^ r ) (\4>o) nas the same phonon occupancies as 
\<j>o) and is involved in the ZPL component of the transition), 
and uj is the frequency of the g-eo ZPL optical transition. The 
energy origin is from now taken at the |eo, 4>o) level. 

The remaining task is the calculation of Pzif). The diago- 
nal electron-phonon matrix elements after the Huangh-Rhys 

transformtion reads A/ e ° e ° = and M e ^ = M eiei - 
M e ° e °. As the calculation of Pz(t) will be made in the new 
phonon mode basis, \<po) will be denoted \4>q) in the following 
in order to keep the notations simple. For t > 0, the evolution 
operator can be expressed as: 



iHt/h _ 



1 

2l7T 



-4- 00 



dEe 



-iE- 



~' h G{E) 



(7) 



where G(E) = 1/ (E + — H) denotes the retarded Green func- 
tion (E + = E + it] with 7] — > + ). In order t o evalu ate this 
term, we use the method of projection operatoiG 5 ! 16 ! 17 ! Let us 
recall a property of this method. If P is a projection operator 
on a subspace of interest and Q = 1 — P, the restriction of the 
Green function G to the subspace defined by P reads: 



PG{E)P 



P 



PH P - PR{E)P 



where the self energy operator R is given by: 

Q 



R{E) = V + V 



E+ - QH Q - QVQ 



V 



(8) 



(9) 



This property will be used twice in the following. First, P 
is defined as the projection operator on the \(p ) phonon state, 
while Q is defined as Q = 1 — P. As the energy origin is 
taken at the |eo, 0o) level, we have: 



(e O) 0o|G(.E)|eo,0o) 



1 



E — D(E) + il(E) 



(10) 



where 1(E) — 3m(eo, (j>o\R(E)\eo, cf>o) (3m stands for imag- 
inary part) and D(E) = ±V f dE'I(E')/(E - E'). The 
energy displacement D(E) is usually found to be small for 
acoustic phonon-electron interaction in QDs and can be ne- 
glected compared to the energy separation between the QD 
electronic levels that are considered, so that D(E) = will 
be taken in the following. Also note that the non-diagonal 
coupling (eo<fio\R\ei(f>o) is found to vanish in the following. 
In addition, we will have often dl/dE < 1 so that 1(E) can 
be replaced by 1(0) = TpinEq. 10 i.e. the ZPL lineshape can 
be approximated by a Lorentzian. Hence the polarization cor- 
responding to the ZPL will exhibit a mono-exponential decay 
of the form P z (t) = e - Tat ' h . 



A. One-phonon processes. 

If we consider only one phonon processes, the term QVQ 
cancels in Eq. [9] so that 



lod(E) 



r = ~f od (-A) 

{(N(E) + l)T od (E) if£>0 
\N(-E)T od (-E) if£<0 



T od (E) 



\M, 



e ei |2 



\ 2 S(E- £l ) 



(11) 



(12) 



(13) 



where N(e) — l/(e 6 ' kBT — 1) is the Bose occupation 
number. The decay rate To/h of the ZPL polarization cor- 
responds here to the Fermi golden rule for acoustic phonon 
absorption (Tq = N(A)T od (A) if A > 0) or phonon emis- 
sion (T = (JV(-A) + l)r od (-A) if A < 0). The population 
decays with the rate 2To / h. 
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B. Two-phonon processes. 

In order to take into account both real and virtual 
transitions, we have to treat simultaneously one and two 
phonon processes. In the calculation of Pz(t) = 
(eo, </>o|e -lff */ fi |eo, (f>o), we will restrict the Hamiltonian to 
phonon states that differ by 2 phonons or less from | <^>o) • 
Within this framework, the calculation will be made non- 
perturbatively. This is achieved in the following using dif- 
ferent projectors. In order to develop Eq.[9] we define a new 
retarded Green function 



G'{E) = 



Q 



E+ - QH Q - QVQ 
Introducing two completeness relation in Eq.|9]gives: 

(eo, MR(E)\e , cj) )=J2Yl V 0a VwG' ab (E) 

4>a 4>b 



(14) 



(15) 



where V a0 = (e , (p a \V\ei, ^p) and G' aj3 (z) = 
{e\,(j) a \G'{z)\ei,(j)p). The term Vba (Ho) is non-zero only 
if (j> a (4>b) differs from <j>o by only one phonon (absorption or 
emission). Hence each 4> a (<f> b ) contributing to the summa- 
tion is either of the form | + i) = aj~ \<fr ) / y ((j) \aiaf\<j)o) or 

| — i) = Oi|0o) / \J ('t'oWt a i\4>o)- This summation involves 
diagonal as well as off-diagonal terms for G', which can be 
expressed separately: 



( eo > |i?(£)|e >o) = J2\V 0a \ 2 G' ai 



^ ^ V a V b *G' ab 
(16) 

We have now to calculate the matrix elements (diagonal and 
off-diagonal) of the operator G' = l/(E + - QH Q- QVQ). 
This operator is similar to the original Green operator except 
that Hq and V have been replaced by QHqQ and QVQ. Sim- 
ilarly, we define new projection operators: P' is the projection 
operator on the subspace defined by by phonon states differ- 
ing from |0o) by one phonon ({</>,; such as Vbi ^ 0}) and 
Q' = 1 — P' . Applying Eqs.[8]and[9]once again leads to: 



P'G'(E)P' = 



F 

E+ - P>H a P> - P'R'(E)P' 



(17) 



at least 2 phonon from the initial state \4>q). As the calculation 
is made up to 2 phonon processes, Q" = P" where P" is 
the projector on states differing by exactly 2 phonons from 
\4>q) (states differing by 3 phonons or more from |</>o) are not 
taken into account). As a consequence the term Q"VQ" in 
G" vanishes since V couples states that differ by one phonon. 
The diagonal matrix elements of R' read: 



E 

i/) 



VaffG'j 



(19) 



where V af = (e u <f> a \V\f) and G" f = (f\G"\f). The |/) 
states differs necessarly by 2 phonons from \4>o). The non- 
diagonal matrix elements of R' read: 



R'nl 



I/) 



VafV fb G" f 



(20) 



In order V a fVf b to be non-zero in this equation, the phonon 
state of |/) should be of the form | ± i,±j) given that <j> a 
and cj) b are of the form | ± i) and | ± j) (|/) = |eo/i, 4>o) is 
forbidden by the Q" projector). 

The non-diagonal matrix elements of R' are of the order 
of 0(1/N cr ) since they involve a product of two terms in V, 
while the diagonal matrix elements of R' involve a summa- 
tion over the N cr phonon modes and are consequently of the 
order of 0(1). In the following, we decompose P'R'P' into 
diagonal R' d and non-diagonal R' od contributions. In the ma- 
trix elements of G', we need to keep only terms up to the first 
order in R' od (i.e. in 0(\/N cr )) since the higher order terms 
give a vanishing contribution in Eq.[T6] Hence inverting equa- 



tion 17 gives: 



P'G'P' = Gl, 



Kj d sx od Kj d 



(21) 



where G' d = P'/(z - P'HqP' - R' d ) is the diagonal part of 
P'G'P'. The non-diagonal elements of G' read: 



f D> fit 



(22) 



where R' a g = (ei, (j> a \R'\ei, tfip) Following the assumption 
that energy displacements are negligible, we keep only the 
imaginary part of G'j(E) ~ —m8(E — Ef). We keep also 
only the imaginary part in the matrix elements of R'. The 
imaginary part of G' aa can then be expressed as: 



P'R'(E)P' = P'VG"{E)VP' 



(18) 



where G"(E) = Q"/(E+ - Q"H Q" - Q"VQ") with Q" = 
}' = 1 — P — P'. Q" is the projector on states that differ by 



3m (G' aa ) 



\G' aa \ 2 Y,\ v *f\ 2G "f ( 23 > 



Now we can express the imaginary part of the self-energy: 



ir 



2 c; 



4>o. <t>b^a f 



V a fVf b G' bb V b Q ) G"f 



(24) 
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where Die stands for real part. The phonon part of | /) dif- 
fers from the initial phonon state |</>q) by 2 phonon occupan- 
cies. |/) being of the form |e /i,±i,±j) (contributions of 
the form |e /i , ±2z) gives vanishing contribution in the limit 
N cr — * oo), one can always define the two state \<f> a ) = | ± i) 
and \cj) b ) = | ± j) so that \e\,4> a ) and |ei,0&) are the only 
states coupled to |eo,<^o) an d |/) simultaneously. As a con- 
sequence, the double summation over the phonon modes \<j) a ) 
and | <t>b) can be replaced by a simple summation over the |/) 
states: 



T = nJ2 \VoaG' aa V af + V Qb G' b * b V bf \ 2 6(s f ) (25) 
/ 

The above formula can be expressed as a sum of two con- 
tribution,which corresponds respectively to the electronic part 
of |/) being either e or ex'. 



1 -1 + L o 



(26) 



We define V° d = (e ,(f>o\V\ex,(l> a ) and V d = 
(ei , 4>o\V\ei,(f> a ). If |/) is of the form |eo, ±i, ±j), we have 
V 0a = V° d , V af = V b nd \ V ob = V b od and V bf = V^ d \ so 
that T ( a) reads: 



r o Q) = lT,T,\ V a 0d \>b°r\GL + Gi* b \ 2 S(e a + e b ) (27) 

The Iq term corresponding to |/) in the e\ electronic state 
reads: 



|ei,#o) 



|eo, 0o) 




eo,-5 +;/) 



FIG. 1: (Color online) Schematic of the electron-phonon levels in- 
volved in the dephasing processes. Levels are classified in three 
columns according to the number of phonon occupancies they differ 
from 4>o' the projector operator P, P', and P" corresponds respec- 
tively to zero, one and two-phonon difference from the initial phonon 
state. The levels in solid lines illustrate coherent couplings involved 
in two-phonon virtual processes. 



where T od (E) was defined in Eq. 13 and 



T d (E) = it Wt x 61 - M^\H(E - Ei) (33) 



^ ) = \ E E ^a^G'aa + V b ° d V?G& | ' S( A + S a +S b ) 

(28) 

To evaluate these terms we need first to express the diago- 
nal matrix elements of R' which are involved in the diagonal 
matrix elements of G' . They can be expressed as the sum of 
two contributions: 

Ka(M + ) = -Hod(sa) - *7d(-A - e„) (29) 
where the following quantity are defined: 

lod (E)=irY,\V a ° d \ 2 6(E~e a ) (30) 



The diagonal matrix elements of G' involved in Eqs. 27 and 
28 can now be expressed as: 



G " a(0) -A - e a + i lod (~e a ) + i ld (~A ~ e a ) (34) 



In the following, we will consider first only off-diagonal 
interactions, and then both diagonal and off-diagonal electron- 
phonon interactions. 



II. OFF-DIAGONAL INTERACTIONS 



ME) 



nJ2\V a d \ 2 S(E-s a ) 



(31) In this section, we consider only off-diagonal interaction in 
V e . p h- In this case, only the term arising from remains: 



7o<2 and j d corresponds respectively to the broadening 
due to non-diagonal and diagonal interaction within the ei 
state. Using the Bose-Einstein distribution function N(E) = 
1/ {e E l kT — l), they can be expressed as 



j(N(E) + l)T od/d (E) if£>0 



lod/d{E) \N{-E)T od/d \-E) 



if E < 



r ° = \ E E S ^ + ^\ V a° d \ 2 \V b od \ 2 \G' a + G b \ 2 (35) 

a b 



(32) Using Sb — — e a , we obtain 
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r °~2 ( ^ ° ((A + e Q )^(7 ,(-ea)) 2 )((A- £a ) 2 + (7o£i ( £Q )) 2 ) (36) 

In the continuum limit for phonons (N cr — ► oo), inserting an integral over the energy leads to: 

r 1 ^ <- ^ 4A 2 + ( 7orf (- g )- 7orf ( £ )) 2 

In order to explicit the temperature dependence, j d can be replaced by its expression as function of T d and the Bose 
distribution function N(e). 

+°° a A 2 + r 2 OV4 

deN(e)(N(e) + l)r 2 d ( £ ) T ; ;/ ^ (38) 



o 



(A + e)3 + [iV( £ )r od (e)] 2 j [(A - e) 2 + [(iV(e) + l)r od (e)] 
As A is usually large compared to T d(e), a simplified expression can often be used (e.g. in ref.fTTTl: 



To = - r + , (39, 

* Jo (A + E) 2 (A - e)> + [(JV(c) + l)r„i(e)] 2 



Physically, this expression contains both the effect of real transitions due to one phonon processes and virtual processes 
involving two phonons. As discussed below, the peak around e — A will be attributed to one-phonon absorption process, while 
the contribution of non-resonant phonon (e away from A) will be interpreted as virtual transitions, i.e. an elastic diffusion of 
phonons without change in the electronic state. The term 4A 2 /(A + e) 2 appearing in the above equation can be understood as 
an interference term between two different path contributing to virtual processes: the intermediate state can be obtained either 
from phonon absorption or emission (see Fig. [TJ. These two path interfere constructively (destructively) for phonon energies e 
smaller (larger) than A. 

III. DIAGONAL AND OFF-DIAGONAL INTERACTIONS 

Including the diagonal interactions, the Tq term is given by: 

r («) 1 [ + °°^ M ( A 4A 2 + ( 7m ( £ )- 7m (-e)) 2 

where 7m (e) = 7 od(£) + ld{— A + e). The term 7m (£) plays a significant role only for e ~ ±A, which involves 7 <j evaluated 

around where it vanishes. As a consequence is almost not affected by the inclusion of the diagonal interactions. 

The term T { P) (which was vanishing in the absence of diagonal interactions) corresponds to real transitions assisted by 2- 
phonon absorption/emission, and depends on the quantity: 

\V a ° d V b d G' a + V b od V d G' b *\ 2 = 2J2J2 \Va° d \ 2 \V b d \ 2 \G' a \ 2 + 2 ]T £ me(V° d V b d G' a V b nd *V d *G' b ) (41) 

a b a b 

In many cases, such as for fundamental intraband or interband transitions in QDs, the second term vanishes because of the 
symmetry of the wavefunctions J2 q M^ l£l M^° ei * = 0. In such case, Ig reads: 



Inserting integrals in the above expression leads to: 



r W> _ I [ + °° dE lo d (sh d (-A-s) 



TV 



do 
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FIG. 2: (Color online) Full width at half maximum (FWHM) 2r 
(full line) of the ZPL component of the s-p x intersublevel transition 
as a function of the energy detuning A pp between the p-levels for 
T=100K. The Fermi golden rule for phonon absorption 2y d(— A) is 
shown in dotted line. The contribution of real 2-phonon transitions 
2Fq appears in dashed line. 



Note that this formula can be written as / c (— A), where 
f c (E) = J dej d(£)A s (E—e) is the convolution between 7 G[ ; 
(Fermi golden rule for one-phonon absorption/emission) and 
the function ^4 s (e) ~ 7d(e)/e 2 describing the one-phonon 
phonon contribution to the sidebands spectrum of the erj-ei 
transition. In analogy with acoustic phonon sidebands of op- 
tical transition, the above formula can be interpreted as the 
apparition of phonon sidebands in a acoustic phonon assisted 
transition. As our calculation is made up to two-phonon pro- 
cesses, the sidebands taken into account involve only one- 
phonon process. For optical transition, however, the appari- 
tion of sidebands does not change the radiative lifetime since 
the relative variation of photon energy is negligible on the 
energy scale of acoustic phonons. On the contrary, acoustic 
phonon sidebands modify the efficiency of acoustic phonon 
induced transitions through an additional contribution of mul- 
tiphonon processes. 



IV. APPLICATION: DEPHASING OF INTERSUBLEVEL 
TRANSITION 

In order to illustrate the theory presented above, we study 
the dephasing of intersublevel transition in InAs/GaAs self- 
assembled quantum dots between the ground s-like state and 
the p-like first excited state (denoted p x here). The higher 
energy p-state is denoted p y . The contribution of acoustic - 
phonon to the dephasing is calculated according to the Eqs. 40 
and [42] In Ref. \TT\ the temperature dependence of the 
linewidth was presented, taking into account 2-phonon pro- 



cesses induced by off-diagonal interactions (Eq. 39 1. Here we 



FIG. 3: (Color online) Linewidth (FWHM) 2r of the ZPL of the 
s-p x intersublevel transition versus temperature for different values 
of the energy splitting A pp between the p-states (the values of A pp 
are indicated on the figure). 



In fig. [2] we present the linewidth of the s-p x intersublevel 
transition as a function of the energy detuning A pp between 
the p-levels. The linewidth given by the Fermi golden rule for 
phonon absorption is also shown for comparison. There is a 
striking difference between the two, which demonstrates the 
importance of two-phonon processes. Application of Fermi 
golden rule results in a overestimate of the linewidth in the 
vicinity of its maximum efficiency, and a strong underestimate 
away from this maximum. 

Note that we have checked that the assumption of an expo- 
nential decay of the ZPL is a good approximation in the full 
calculation: at 100K, dl/dE < 0.3. In addition, at the same 
temperature, the energy displacement D(E) remains smaller 
than 0.15 meV. 

In fig. [3] the s-p x ZPL linewidth as a function of tem- 
perature is plotted for different values of the p x -p y splitting 
App. It is worth to mention that the present calculation does 
not include the populati on relaxa tion from p x to s by an- 
harmonic polaron dec ay^^- 8 l 1 9 I 2 ° I 2 1 1 which contribute addition- 
ally to the linewidth, and is the dominant mechanism at low 
temperature 11 . Note also that the ZPL component in the in- 
tersublevel absorption decrease with increasing temperature 
down to approximatively 0.5 at room temperature. 

A physical interpretation is provided by studying the inte- 

CTg = / deh(e)). h{e)de cor- 
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take into account additionally the contribution of 2-phonon 
real transitions due to diagonal terms (Eq. [42}. The QD pa- 
rameters are the same as in Ref. Qj] with a s-p x transition 
energy of 53 meV. 



grand (denoted h) in Eq. 
responds to the contribution" to the dephasing (non-diagonal 
terms only) due to phonons with energy comprised between 
e — de/2 and e + de/2. In Fig. [4| the function h is plotted for 
different value of the p-p splitting. The function h is sharply 
peaked in e = A, which corresponds to real transition involv- 
ing one-phonon process, i.e. phonon absorption promoting the 
electronic state from e to e.\. If the condition d^^jdE <C 1 
is fulfilled (cases c and d), the area of this peak can be approx- 
imated by j d{— A) which corresponds to the Fermi golden 
rule for phonon absorption. In addition, in the cases a, c and 
d of Fig. [4j the function h presents an additional broad peak 
between 1 and 4 meV which is attributed to 2-phonon virtual 
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made non-pertubatively within this framework. An analytic 
formula has been given in terms of the electron-phonon ma- 
trix elements. The theory has been illustrated to the study of 
the dephasing of intersublevel transitions in QDs. The con- 
tribution of phonons to the dephasing has been studied as a 
function of their energy. 
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APPENDIX A: ACOUSTIC PHONON SIDEBANDS 



FIG. 4: (Color online) Function 2h representing the contribution 
to the linewidth 2Tq as a function of the acoustic phonon energy 
(T=100K) for different p-splittings. (a) lmeV (b) 3 meV (c) 5 meV 
(d) 7 meV. The area under the curve (red) is equal to 2Fq (FWHM). 



processes. 

In fig. [4]b, the contribution to one and two-phonon pro- 
cesses cannot be separated. In this case, the calculated 
linewidth is found to be smaller to the Fermi golden rule (see 
Fig. |2j. This feature can be explained qualitatively by replac- 
ing delta energy conservation S(A — e) appearing in Fermi 
golden rule by Lorentzian of linewidth 7 0( ;(A). As a conse- 
quence the strong variations in j d(— e) are smoothed (over 
an energy 7 d(£)), resulting in a negative correction for the 
peak of j d(-e). 

Depending on A pp , the relative contribution of real and vir- 
tual transitions varies. In Fig.|4]c, for A pp — 5 meV, the area 
of the two peaks are of the same order. For larger p x -p y split- 
ting (Fig. [4}i), dephasing processes are dominated by virtual 
transitions. This can be explained by noticing that the phonons 
which are the more efficiently coupled to the electronic states 
have wavelength of about the dot size, which corresponds to 
typical energy of 1 or 2 meV in self-assembled QDs. This can 
be seen by looking at the variation of 7od( _ e) (Fermi golden 
rule for absorption) on fig. 2. In case phonons with QD wave- 
length are not resonant with the electronic transition energy, 
they cannot induce real transitions, but can still be elastically 
scattered by the QD and consequently induce decoherence. 
On the contrary, the efficiency of real absorption processes 
decreases rapidly when increasing A pp . This qualitatively ex- 
plains why virtual transitions can be the dominant dephasing 
mechanism. 



Conclusion 

We have presented a non-pertubative calculation of acous- 
tic phonon induced dephasing of optical transitions in QDs. 
The Hilbert space has been limited to states differing by up to 
2-phonons from the initial state, and the calculation has been 



In order to calculate the polarization decay, we decompose 
the initial phonon state \<fio) on the new phonon eigenstate 



(e 0) ^ \e- tHt ^\e , 0o> - ]T |(0 o |^)| 2 (e o , <Pi\e- im ^\e , </>. 

iHt/h 



^2 S (^o|&)(<M<ft))( e o,<fol e ~ ">■"// 



(Al) 



The second term is of the order of 0(1/N cr ) and thus van- 
ishes. The states \4>i) involved in the above equation differs 
from the state |0o) (corresponding to the zero-phonon line) 
by a finite number of phonon occupancies, so that in the limit 
N rr — > oo we have: 



(e ,^\e- im ^\e ,$i) 



<eo,0o|e-^ t/fi |eo,<M 



This allows us to write: 

(e ,Me- lHt/h \eoAo)=Pz(t)g(t) 

P z (t) = {e ,Me- iHt/H \e ,To) 



eo)t/H 



(A2) 
(A3) 
(A4) 
(A5) 



The function g{t) defined here is the Fourier transform of 

g{E) = \(<t>o\$i)\ 2 8(E - (e~ - e )), (A6) 



which is the function describing the acoustic phonon side- 
bands absorption as a function of the energy detuning with 
the zero-phonon line. ^^^^ 
As result from the independent boson model theory^l22l2il 
g(t) can be expressed as: 



g(t) = exp 



/+oo 
dEf(E) ( e - lEt/h - 1 
-OO 



(A7a) 
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(A7b) 

\(N(E) + l) if£>0 
X [N(-E) if£<0 

The function g(t) describes the evolution of the coherence due 



to phonon sidebands only, i.e. taking into account diagonal 
electron-phonon interaction only. The acoustic phonon side- 
bands have been studied extensively in the literature for inter- 
band transitions in both energy domain 4 5 and time domain 6 , 
and more recently for intraband transitions^. 
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